Generation of spin currents via Raman scattering 



in 
o 
o 

wo: 
< 

On 



i 

C 

o 
o 



> 
o 
m 

(N 
00 

o 
in 
o 
-i— > 



i 

C 

O 

o 




AZ=[11 



,z=[110] 

-►jHHO] 



Ali Najmaie/'B E. Ya. Sherman, 1 and J. E. Sipe 1 

1 Department of Physics and Institute for Optical Sciences, 
University of Toronto, 60 St. George St., Toronto, Ontario, Canada M5S 1A7 

(Dated: February 2, 2008) 

We show theoretically that stimulated spin flip Raman scattering can be used to inject spin 
currents in doped semiconductors with spin split bands. A pure spin current, where oppositely 
oriented spins move in opposite directions, can be injected in zincblende crystals and structures. 
The calculated spin current should be detectable by pump-probe optical spectroscopy and anomalous 
Hall effect measurement. 
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Spintronics is the strategy of using the spins of carri- 
ers in solid state structures as a new degree of freedom 
to carry and transform information 0. An important 
capability is the generation of spin currents, which cou- 
ple spin and charge. Suggestions for the generation of 
spin currents have included proposals based on the use 
of the extrinsic and intrinsic spin Hall effects, as well as 
spin pumping [3- 0.0. All-optical generation and control 
of spin currents have also been proposed and observed 
H, IE SHI El- I n these optical schemes, the generation of 
the spin current is accompanied by the deposition of pho- 
ton energies of about 1 eV per carrier in the system. Most 
of this energy is used simply to generate the electron-hole 
pairs, and is lost in the sense that it does not contribute 
to the motion of the carriers. An obvious goal for the all- 
optical manipulation of spin is the minimization of such 
energy loss. 

Here we propose a new scheme for generating spin cur- 
rents in doped semiconductors via a two photon Raman 
scattering process 0, ^} , which can be understood as 
the absorption of a photon of frequency w\ and the emis- 
sion of a photon at lj 2 - This process deposits an energy 
of h£l = h (u>i — u>2) per carrier in the crystal, which for 
GaAs is only about 0.5 meV. In the presence of spin split- 
ting of the bands, stimulated spin-flip Raman scattering 
can produce a pure spin current, where electrons with 
opposite spins move in opposite directions. Hence no net 
charge current is injected. 

Crucial to this effect is the spin-splitting of states, 
which results from the lack of inversion symmetry of the 
crystal. A schematic picture of the conduction states is 
shown in Fig. 1, where Ep denotes the Fermi energ 
For simplicity we take the temperature to be zero 
The proposal in this letter is based on optically induc- 
ing transitions between states d(±k) and u(±k) at ±k 
wavevectors, as depicted in Fig. 1, via Raman scatter- 
ing. As an example, we first consider an optical field 
E = E(ui)e~ Mlt + E(u 2 )e~ M2t + c.c. in the material, 
where c.c. denotes complex conjugation and u>i > u>2 
and then present an experimentally relevant scenario 
of pulse excitation. The interaction of this field with 
the crystal in the dipole approximation is described by 




FIG. 1: (Color online) Schematic picture of the spin-split 
conduction bands. Spin-flip Raman scatterings correspond to 
transitions from d(±k) to u(±k) at ±k. The carriers at u(±k) 
then have opposite spin orientations and hence a pure spin 
current is injected in the system. An example of the generated 
spin current is presented where two pure spin currents with 
orthogonal velocities are shown schematically. The dotted 
lines correspond to the velocities of the carriers, and the thick 
arrows depict their spins. 



H^t = — (e/mc)p • A(t), where e(m) is the charge(mass) 
of an electron, c is the speed of light and A(t) is the vec- 
tor potential associated with the field. In the ground 
state |0), electrons occupy states up to the Fermi level. 
The wavefunction of the perturbed system at time t is 
then written as I \I> (t)) = cq (t) 1 0) + 7 Ck (t) I ehk) , where 

^ — 'k 

|e/ik), denotes the state where an electron (e) and hole 
(h) are generated at point k (see Fig. 1). Time depen- 
dent perturbation theory is used to evaluate Ck(t) and 
hence to describe the rate of change in the expectation 
value of a single particle operator 9 in the independent 
particle approximation, 
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This leads to the evaluation of injection rates of interest, 
including that for the density of flipped spins and for the 
spin current density |l3j . The former can be written as 

^ = e hc \ui,m^i)E b *{u2)E c *{ui)E d {u> 2 ), (2) 
where f &cd (wi, to) = V^ 1 V T abcd {k)S(fl-cj eh (k)) and 
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Here w nm = w„ — £j m , is the normalization volume, 
a,b,c and d are the Cartesian indices and v nm (k) is the 
velocity matrix element between single particle states 
n and m at k. Repeated indices are summed over. 
To quantify the spin current density J a b, we use the 
conventional [Til flsT ] single particle spin current oper- 
ator j ab = (l/2)(v a S b + S b v a ), where S = ficr/2 and 
er = (<r x , a y , cr z ) are the Pauli matrices. The spin cur- 
rent density injection rate is 

_ (4) 
where p abcdf9 {ui, to) = V^V 6j ab T cd f3{k)5(to ~ 

— 4 ' ehk. 

uj eh (k)) and Sj a b = (ek\j ab \ek) - {hk\j a b\hk} is the spin 



current contribution per spin flip [ljj. Both the fourth 
rank tensor £ abcd = £ cdab anc [ sixth rank pseudotensor 

^abcdfg _ ^abfgcd can ^ e s J 10wri £o be rea J by us i n g time 

reversal symmetry; the Td symmetry group of the crystal 
can be used to deduce the non-zero and independent ten- 
sor components. To calculate the tensors, we describe the 
spin splitting of the conduction states by the Dresselhaus 
Hamiltonian in the T point basis { 1 1/2, 1/2) , 1 1/2, - 1/2) } 
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H = 



2m r 



oldK ■ er, 



(5) 



where k x = k x (ky — k%) and K VtZ are cyclic permuta- 
tion of n x , up, — 27.6 eVA 3 is the Dresselhaus param- 
eter and m c = 0.067m is the electron effective mass 
[l?! lrij |. The spin splitting at wavevector k is 2aDH. 
The valence states of the semiconductor serve as the in- 
termediate virtual states n in Eq. ©. The heavy and 
light valence states are modelled using the isotropic 4x4 
Luttingcr-Kohn Hamiltonian, with the Luttinger param- 
eters 7i = 6.85 and 7 = (73 + 73) /2 = 2.5. The split- 
off band is described by a diagonal 2x2 Hamiltonian, 
with the effective mass taken to be 0.154m [J^- A band 
gap energy of E g = 1.51 eV and a split-off energy of 
0.34 eV are used. The Fermi energy of the electrons in 



the conduction band is given by Ep = Eg + Ti 2 k F /2m c , 
where kp = (S^N) 1 / 3 for the carrier concentration N; 
we use a concentration of N = 10 18 cm -3 in the cal- 
culations. The speed of the injected carriers is given 
by vp — hkp/m c « 500 km/s. The intraband velocity 
matrix elements are evaluated using v = Ti~ l dH (k)/9k 
as the velocity operator. We use a Kane parameter 
P = (h/m)(S\ Px \X) = 10.49 eVA and numerically 
evaluate the tensor components. 

Spin orbit coupling splits the Fermi-surface into two 
surfaces. The fraction of the electrons that are con- 
tained in the region between the two surfaces is given 
by w 2m c akp/h 2 and is about 1% of the free carriers 
available in the system at TV = 10 18 cm~ 3 . Only this 
fraction of the free carriers in the system is subject to 
Raman scattering. This sets a limit of « 10~ 2 iV on the 
concentration of carriers that can contribute to the in- 
jected current. 

Spin flip Raman scattering has typically been studied 
using orthogonally polarized fields E^i^) 01- Here- 
after we use the following definition of the coordinates, 
x = [001], y = [110] and z = [110] and consider more 
general fields normally incident on a crystal with a [110] 
growth axis 



E(wi, 2 ) = £a,2/3i,2 (cos(pi 2 )^ + i cos(p? 2 )y) 



(6) 



for phases (pf, p\, pf, pf), and normalization factors 0\$. 
This allows the most general analysis of the fields polar- 
izations. We find the spin-flip and spin current density 
injection rates 
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where g = (0x02) (cos(pf ) cos(pf ) + cos(pf ) cos(p|)) , 
and 7i t 2 are the incident intensity of the £i ;2 = Tv^i,2 
components of the fields. Here A ei , B £l and C ei are linear 
combination of the tensor components £ abcd and ^ abcd f9. 
The mismatch in the dielectric function e across the in- 
terface leads to a ratio 2/(y/e + l) of the electric fields in 
the crystal and air. Two pure spin currents are injected 
in the crystal, one injected along y = [110] with spins 
parallel to z = [110] and another injected along z = [110] 
with spins pointing parallel to y = [110] (Fig. 1). Clearly 
for either right (or left) circularly (e.g. 0i= 02 = 1/V2, 
Pi = P2 = 0' Pi = p\ = 0( or 7r )) an d cross-linearly (e.g. 
(h = 02 = 1, Pi = p\ = 0, pf = pf = 7r/2) polarized 
fields, the same spin currents are injected. Moreover, for 
co-linearly polarized fields, where g = 0, no spin-flips and 
hence no pure spin current can be induced. 

First we consider the case where the photon energy 
Ei = fiwi = 0.754 eV is close E g /2. All two photon 
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FIG. 2: The thin lines show the coefficients, Ai, B\ and Ci 
Eq. |0 for Q. = uj\ — u)2 and Tiuj\ — 0.754 eV, about half the 
band gap. The thick lines show the coefficients, An, Bn and 
Cn Eq. Q for JkJi = 1.49 eV, about 20 meV below the band 
gap. 



absorption from the valence to the conduction bands is 
avoided in this scenario. In Fig. 2, we present A\, B\ and 
C\. The role of the resonance enhancement can be inves- 
tigated by using optical fields of energy En = hui 1 = 1.49 
eV close to the band gap. Transitions from the partially 
filled conduction band to the higher bands are avoided 
at this energy, but two photon excitations from the va- 
lence states to conduction states are allowed. Nonetheless 
these transitions are away from the T point and should 
not contribute to the injected pure spin current from the 
spin flip Raman processes we discuss in this letter. The 
two photon absorption is estimated to be an order of mag- 
nitude smaller than the spin-flip injection rate discussed 
here [2fj , and excitons are not generated. The carrier and 
spin current density injection rates characterized by An, 
B\\ and Cjj are presented in Fig. 2. The sharp features in 
the results, shown in Fig. 2, correspond to the van Hove 
singularities in the joint density of states for the spin-flip 
transitions. The van Hove singularity is due to the van- 
ishing derivative in the spin splitting, 2cv£>VkK along high 
symmetry directions, and an analysis of the Dresselhaus 
Hamiltonian shows that the resulting density of states 
increases as ln(l/ |f2 — fl V H\), where Q v h is the position 
of the singularity. To this point we have discussed the 
injection rate of pure spin currents by stimulated Raman 
scattering using two nominally continuous wave (CW) 
beams at distinct frequencies u>2 and loi > lu2- We now 
turn to the ultrafast pulse regime, and stress that a pure 
spin current can be injected through the Raman effect 



by a single pulse, with the scattering now arising from 
the spectral width of the pulse. For a pulse with a centre 
frequency loq and width Aw, if we have ujq 3> Alu and 
luq fimax, where 7tf2 ma x is the maximum spin splitting, 
the straightforward generalization of the calculation in 
Eqs. J3J) and (0} leads to the injection of densities of spin 
flips and spin current given by 
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7 g r ^( W0; tt)G afccrf (f>), 
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where the field correlation function for E(i) 
(27r)- 1 /_^ o dwE(w)e- l " t is 



G abcd {VL) 

OO 

duJiduj^ 



(10) 



E a (u 1 )E*>*(u 1 - Q)E c *(u' 1 )E d (Lj' 1 - O). 



We consider a (right or left) circularly polarized short 
pulse of the form E(t) = 2E a exp(— t 2 /r 2 ) cos(woi), with 
fiuence (pulse energy per unit area) F. A short pulse [2l| 
max 1 and r <C (iTn c /\e\, where u is the mobility) 
essentially leads to the generation of the spin current in 
the ballistic limit. For hujQ — 0.754 eV (corresponding 
to Aj, B\ and C\ in Fig. 2), the density of spin flips 
per pulse is Ani ~ (1.8 x 10 20 cm/ J 2 )F 2 , which gives 
Anj w 10 15 cm -3 for an experimentally accessible fiuence 
of_F = 3xlO~ 3 J/cm 2 . In a scenario where Tiujo = 1.49 
eV (corresponding to An, B\\ and Cn in Fig. 2) the spin 
flips density due to Raman scattering is Ann ~ (1-2 x 
10 22 cm/J 2 )F 2 ; a spin flip density of Ann w 10 15 cm~ 3 
can be achieved for a fiuence of F = 3 x 10~ 4 J/cm 2 , 
which reflects a resonance enhancement of the process 
by two orders of magnitude. 

Note that the injection in this process is proportional 
to F 2 , whereas from Eq. Q for nominal CW irradiation 
one might expect a scaling with / max r cx F 2 /t where 
J max is the maximum intensity in the pulse. The reason 
is that, for O max r <C 1, a factor r2 max T w f2 max /Ati; 
appears in the final result following from Eqs. JSJ) and 10, 
reflecting the fact that only a fraction f2 max /AcL> of the 
pulse bandwidth can contribute to the Raman scattering. 

With each pulse a number of carriers ~ An x S-q are 
excited, where S is the laser spot area and n is the pen- 
etration depth. In earlier experiments 0, |g involving 
interband transitions, a large number of carriers are avail- 
able in the valence states, and typically ~ (10 18 cm _3 )S'n, 
where r\ w 1 /mi, carriers can contribute to the spin cur- 
rents. On the other hand, in Raman injected spin cur- 
rents the penetration depth is much larger than that for 
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interband transitions, with a limit rj « 10 3 y«m set by 
the Drude absorption due to free carriers. This allows a 
large ss (10 15 cm _3 )5'?7' number of spin flips in the sam- 
ple, which will then contribute to the injected pure spin 
current. 

Pump and probe schemes have been used to exper- 
imentally observe injected pure spin currents 0. An 
experimentally relevant parameter in these schemes is 
the separation between the opposite spins after momen- 
tum relaxation. This distance can be estimated by 



H.M. van Driel, and A.L. Smirl for useful discussions. 



d a = (2 M m c /| e |n)E„ 



J ab 



I 



For circularly 



polarized as well as cross linearly polarized fields, a spin 
separation of d v « 18 nm for fi w 5 x 10 3 cm 2 /(Vs) is 
predicted, a distance which is within experimental ca- 
pabilities 0. Measuring a pure spin current by using 
the anomalous spin Hall effect arising due to skew scat- 
tering by impurities 0> Ell nas also been proposed || . 
The spin current J a b causes a spin-Hall bias V s h along 
the direction perpendicular to the electron propagation. 
This bias is estimated as V s h ~ tan(# s #)14ff, where 9 s h 
is the spin-Hall angle and V e s is the effective bias that 
would cause a current density eJ a b/h', it is of the order 
of eAnvF, where An is the number of carriers that form 
the pure spin current. Therefore V^ff ~ L (An/N) vp/fj,, 
where L f» 10 /jm is the lateral size of the system, which 
is assumed to be approximately equal to the spot size. 
The ratio An/N is estimated as (aokp/ Ep, where £ 
is the fraction of Raman accessible electrons that un- 
dergo a spin flip. With the above mobility we obtain: 
V e s/L w 10( V/cm. The above analysis shows £ to 
be about 0.1. With the angle 8 s h estimated by Abaku- 
mov and Yassievich [24| for GaAs at low temperatures as 
9 s h ~ 10~ 3 , we obtain an upper estimate of V s h ~ 10~ 6 
V. 

In conclusion, stimulated spin-flip Raman scattering 
can be used to inject spin currents in doped semiconduc- 
tors with spin splitting. We have considered the case of 
an n-doped GaAs bulk crystal, and showed that the in- 
jected spin current is a pure spin current comprised of 
spins with opposite orientations travelling in opposite di- 
rections. This pure spin current should be experimentally 
accessible, and can be achieved either with two contin- 
uous waves contributing to the Raman scattering, or by 
using a single pulse. The scheme is all-optical but, unlike 
other optical schemes involving the excitation of electron- 
hole pairs, here there is minimal energy deposited in the 
crystal. The system is thus only slightly perturbed from 
equilibrium, and the description of transport and the de- 
sign of device structures should be easier than in sce- 
narios employing other optical schemes. High mobility 
quantum wells should provide ideal testing grounds for 
the general proposal discussed in this letter. 
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